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SUMMARY 

This ASTM Round Robin was conducted to evaluate the state of the art 
in stress analysis of adhesively bonded joint specimens. Specifically, the 
participants were asked to calculate the strain-energy-release rate for two 
different geometry cracked lap shear (CLS) specimens at four different 
debond lengths. The various analytical techniques consisted of 2- and 
3-dimensional finite element analysis, beam theory, plate theory, and a 
combination of beam theory and finite element analysis. The results were 
examined in terms of the total strain-energy-release rate and the mode I to 
mode II ratio as a function of debond length for each specimen geometry. 
These results basically clustered into two groups: geometric linear or 

geometric nonlinear analysis. The geometric nonlinear analysis is required 
to properly analyze the CLS specimens. The 3-D finite element analysis 
gave indications of edge closure plus some mode III loading. Each partici- 
pant described their analytical technique and results. Nine laboratories 
participated. 
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INTRODUCTION 


Many applications of adhesives result in bonded joint 
geometries that have complex stress states at the bondline 
termination. Debond initiation and growth may occur which could 
lead to joint failure. Because debonding is a self-similar cracking 
problem, it is natural to describe the phenomenon in terms of 
fracture mechanics. Furthermore, depending on the joint geometry and 
loading conditions, these stress states can result in three modes of 
debond growth: 

(1) Opening mode I, due to tensile stress normal to the plane of the 
debond, results in a Gj strain- energy-release rate. 

(2) Shearing mode II, due to in-plane shear stress, results in a Gjj 
strain-energy-release rate. 

(3) Tearing mode III, due to anti-plane shear stress, results in a 
G III strain-energy-release rate. 

The calculation of the strain energy release rate under such 
conditions is not trivial, especially when geometric non-linearities 
may exist. If standard test methods are to be developed to 
assess mixed mode static fracture toughness and debond propagation 
rates for designing adhesively bonded joints, confidence must be 
obtained in the methods used to calculate the mixed mode values of 
strain energy release rate. By comparing various analytical 


1 



approaches to common problems can one determine which techniques 
give acceptable answers. 

This problem of calculating the strain energy release rate for 
adhesive joints has been address by ASTM Task Group E24.04.09 on 
Crack Growth in Adhesively Bonded Joints. ASTM Committee E-24 on 
Fracture Testing sponsors this task group. The task group was formed 
in 1981 to evaluate and recommend test methods for characterizing 
debond propagation in adhesively-bonded joints. In the beginning 
the task group decided that prior to conducting tests or 
recommending testing procedures it needed to assess the present 
state of the art in stress analyses that can be used to study debond 
propagation under mixed-mode conditions in adhesively-bonded joints. 
For this purpose a round robin was conducted and is described 
herein. 

The cracked-lap-shear (CLS) specimen [1] was chosen for 
evaluation. The CLS specimen has been used for studying composite 
delamination [2, 3, 4, 5] and adhesive joint debonding [1, 6, 7]. 

The debond tip in the CLS specimen is predominatly loaded in mode II 
with approximately thirty percent mode I. The CLS specimen is 
representative of the mode mix found in many applications of bonded 
joints and was therefore chosen as an appropriate specimen for the 
round robin. There is no known "exact" solution for the CLS 
specimen. The round robin results will be compared to each other 
and to experimentally observed specimen behavior to see if any 
consensus "correct" solutions exist. 
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The two CLS specimen configurations analyzed are shown in Figs. 

1 and 2. As shown in the figures, the difference in the two 
configurations is the thickness of the lap adherends. The adherends 
are aluminum. Table I gives the assumed material properties. Boundary 
conditions are also shown in Figs. 1 and 2. A 11.12 kN tensile load 
was applied for both specimen geometries. Perfect bonding between the 
adhesive and adherend was assumed. The material properties were 
considered to be linear elastic with no time-dependent behavior. For 
this round— robin effort, the debond was located at the middle of the 
adhesive layer although that is not usually the case in practice [7]. 
The debonding normally occurs near the strap adherend. 


The information sought from each participant was their 
calculated strain-energy-release rates, Gj and Gjj. The tearing 
mode GUI was reported, if available. Debond lengths of 2.54, 6.35, 
25.40, and 101.60 mm were analyzed for both specimen configurations. 

Table II lists the organizations and people who participated in 
this round robin. Also listed in the table are their techniques 
used to predict the strain-energy-release rates and the symbol 
representing each technique used in subsequent data plots. Table II 
also gives the appendix number for each participant. Each appendix 
gives details about the analysis techniques used and the resulting 
predictions. The techniques fall into four general classes: Closed 
Form, Geometric Linear Finite Element Analysis, Geometric Nonlinear 
Finite Element Analysis, and mixed Closed Form / Finite Element 
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Analysis. All reported analysis used linear elastic material 
properties. In this paper linear or nonlinear analysis refers to 
geometric linear or geometric nonlinear analysis. 

PREDICTED RESULTS 

Figs. 3 and 5 show the predicted values of total strain-energy- 
release rate, Gfj>, versus debond length, for the equal and unequal 
thickness adherend, respectively, where G T = Gj + G I3; + G IXI . The 
plotted values of G^ are as calculated; however, the results are 
grouped according to the analysis method. The first group (going left 
to right) is the Closed Form results: the second group is the 
Geometric Linear Finite Element results; and the third group is the 
Geometric Nonlinear Finite Element results. The third group includes 
the mixed Closed Form/Finite Element Analysis method. A symbol is 
assigned to each prediction technique as given in Table 2. Figs. 4 
and 6 show the predicted ratio of Gj/Gn versus debond length for the 
equal and unequal adherends, respectively. The plotted values of 
Gj/Gjj are grouped as discussed above. 

For the equal thickness adherend specimen results shown in Fig. 

3, the four geometric nonlinear analyses resulted in a rather constant 
G T as a function of debond length, while the geometric linear analyses 
resulted in a "humping" of the data with debond length. There is 
reasonable scatter (less than 10 percent) within each analysis 
approach (geometric linear or nonlinear) . The 3-D geometric linear 
predictions resulted in a G T about 10 percent higher than the 
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average of the 2-D geometric linear. The geometric nonlinear analysis 
results of Law (#) are approximately 10 percent higher than the other 
nonlinear results. (Notice that Law only predicted G T at debond 
lengths of 2.54 and 101.6 mm). Law assumed plane stress conditions 
while the other participants chose plane strain conditions for their 
2-D analysis. The results of Erdogan's Reissner plate theory (A) 
approach are very close to the 2-D geometric linear F.E. analysis 
while Brussat's closed-form beam theory (-)-) method is close to the 
geometric nonlinear F.E. analysis. 

The data trends for the Gp/Gpp ratio predictions of the equal 
thickness adherends shown in Fig. 4 is similar to those previously 
discussed for G T except for the Lof's 3-D predictions 0 ) . Lof's 
G l/ G II ratios are noticeably lower than those predicted by any of the 
other techniques. 

The G<p behavior of the unequal thickness adherend shown in Fig. 

5 is similar to that of the equal thickness adherend. Again the 
plate theory and the geometric linear analyses showed a "humping" 
with increasing crack length while the virtual crack extension/beam 
theory method and the geometric nonlinear analyses gave constant G T 
with increasing crack length. Mall's (®) and Hufferd's 0 analysis 
were about 10 percent below the other two nonlinear analysis and the 
beam theory for some unexplained reason. 

The Gp/Gu ratio predictions shown in Fig. 6 for the unequal 
thickness adherends followed the same trends as the G T predictions 
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with two notable exceptions: the beam theory G x /G „ ratio was much 
higher than that predicted by the geometric nonlinear analyses; and 
the combination beam theory/finite element analyses resulted in a 

significantly higher Gj/Gjj ratio at 101.6 mm than at the shorter 
debonds . 


All of the finite element analyses used Rybicki and Kanninen 
[9] virtual crack closure technique to calculate values of G x and 
G IXl except for Hufferd's and Lof's. They calculated stress 
intensity factor, K, then converted to strain-energy-release rates 

DISCUSSION 


There is excellent agreement among the 2-D geometric nonlinear 
analyses of the equal thickness adherends CLS specimen. The geometric 
nonlinear analysis of Law (*) assumed plane stress conditions while 
the analysis of Everett and Whitcomb ( 0 ) , Mall (®) , and Dattaguru and 
Mangalgiri ((*) assumed plane strain conditions. The plane stress 
analysis will result in G T being 1/(1-V 2 ) (i.e. 1 . 19 ) higher than the 
plane strain analysis. This is approximately the difference observed 
in Fig. 3 . In appendix VII, Law shows how the geometric nonlinear 
analysis requires several iterations to converge. He also discusses 

the importance of modeling the entire bondline and not just a short 
portion. 

For the equal and unequal adherend specimens the beam theory 
predictions of G T are very close to the geometric nonlinear 
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predictions. As seen in Figs. 4 and 6, the values of the geometric 
nonlinear Gj/Gjj ratio is about the same for the equal and unequal 
thickness adherends specimens (0.25 and 0.225, respectively) in 
spite of the rather significant change in geometry. The beam theory 
analysis results were 0.27 and 0.36, respectively. Dattaguru, 

Everett, Whitcomb, and Johnson [8] have suggested that the beam 
theory approach is accurate for Gj/Gjj determination for only equal 
thickness adherends specimens. Brussat did not include the debond 
length, a, in the calculation of Gp, Gp, and Gpp; thus, these values 
are independent of debond length. 

In Appendix I, Anderson, Abrahamson, and DeVries calculate Gp for 
a debond length of 50.8 mm (in addition to the four lengths requested) 
in the equal thickness adherend using a geometric linear analysis. The 
Gp at 50.8 mm was significantly higher than for the 25.4 mm debond 
length. This indicates that the "hump" predicted by geometric linear 
analysis is even higher than shown in Fig. 3, indicating a larger 
discrepancy between geometric linear and nonlinear analyses. In 
Appendix IV, Erdogen and Joseph present Gp as a continuous function of 
debond length, also illustrating the "humping” behavior. 

In Appendix I, Anderson, Abrahamson, and DeVries presented two 
sets of analytical results. The first set used a rather course mesh 
and a Gp of 375 J/m 2 was obtained. The mesh was then refined and a 
value of Gp equal 213 J/m 2 resulted. This was consistent with the 
other 2-D linear finite element analysis. Thus, Anderson, et al. 
found that too course of a mesh could give poor results. 
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The geometric nonlinear analyses presented for the two example 
problems herein predicted G T to be practically constant with debond 
length. This is not always true. For specimens with shorter lengths 
or thicker strap adherends the predicted G T may vary with debond 
length [8,10]. However, the constant values of G T with debond lenght 
are believed to be the correct results for these geometry specimens. 
Debond growth rates have been shown to be constant over crack lenght 
in similar geometry specimens made with adhesives shown to have growth 
rates that are governed by G T [7,11]. Furthermore, in reference 8, 
measured displacements in a cracked lap shear specimen agreed closely 
with the geometric nonlinear analysis and not with the linear 
analysis. All of this points to the need to use a geometric nonlinear 
technique for analyzing the CLS specimen [4,8], 

For those configurations of the CLS specimens that do show a 
variation in the geometric nonlinear G T with debond length [8] , the 
beam theory approach will give incorrect results. This is because 
the beam theory will give only one value of G T since the theory 
assumes an infinite specimen length. 

Mall's (©) approach models the crack tip area with finite 
elements to calculate strain-energy-release rate but uses the beam 
theory for overall specimen behavior. For equal thickness adherends 
specimen. Mall's results are about 5 percent above the geometric 

nonlinear finite element results. For the unequal thickness 

adherend, Mall's results for shorter crack lengths are about the same 


8 



as for the other nonlinear predictions ; however, his results are 
higher for the longer crack length. 

Lof's (Q>) results were very interesting. He used a 3-D 
geometrically linear finite element method. The G<p and Gj/Gjj ratios 
were not very close to the 2-D analyses for reasons unknown. 

However, his analysis showed that Gj and Gjj varied across the width 
of the specimen and that a Gjji existed near the edges of the 
specimen. The Gp remains fairly constant across the width of the 
specimen while the Gj, Gjj, and Gpjj components vary. The edge 
closure discussed by Lof is quite surprising. The effects of edge 
closure, which is obtained only by the 3-D analysis, needs to be 
further studied. 


CONCLUSIONS 

Fracture mechanics is being used to describe static toughness 
and damage growth behavior in adhesively bonded joints. Because the 
stress states at the crack tip in a bonded joint can be complex 
and result in several potential modes of crack propagation, it is 
critical that the joint specimens be analyzed properly in order to 
insure correct interpretation of experimental test results. To aid 
in evaluating potential techniques for analyzing bonded joints, 

ASTM Committee E-24 sponsored a round robin stress analysis of the 
cracked-lap-shear specimen. The participants were asked to 
calculate the strain-energy-release-rate at four different debond 
lengths on two different geometry specimens. The results of this 
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round robin yielded the following conclusions: 


1. There was good lab-to-lab consistency of predictions between 
analytical techniques of the same type. For example, all of the 2-D 
geometric linear techniques gave results within 5 percent of each 
other. These results suggest that although no generally applicable 
closed form solutions exists for the cracked lap shear specimen, good 
consistent results can be obtained by properly using finite element 
techniques. 

2. The geometric nonlinear analysis techniques give results that 
are most consistant with observed experimental behavior and should 
therefore be used for analyzing the cracked lap shear specimen. 

3. Brussat's beam theory approach gives a good closed-form 
approximation of total strain— energy— release rate for the example 
CLS specimens. However, beam theory does a poor job of predicting 
G l/ G II ratio when the adherends are of different thickness. 

4. Care must be taken to ensure proper modeling of the specimen in 
order to get correct results. For example, the adhesive bondline must 
be modeled with proper grid density (as discussed by Anderson, et al.) 
Also, the full bondline, not just the near crack-tip region, needs to 
be included in the model (as discussed by Law) . 
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5. The three dimensional analysis showed that an edge closure 
effect may be present at the debond front of the cracked lap shear 
specimen. This analysis also indicated that a small amount of Gjjj 
was present near the edge. The G«p remains nearly constant across 
the debond front although the mixture of Gj, Gjj, and Gm varied. 
This three-dimensional effect discussed by Lof needs to be studied 
further. 
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Table I. - Material Properties 



Aluminum adherend 

Adhesive 

E 

72450 MPa 

1932 MPa 

V 

0.33 

0.40 
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TABLE II - ROUND ROBIN ANALYSES & PARTICIPANTS 


ANALYSIS 

NAME 

ORGANIZATION 

TECHNIQUE/MODEL 
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APPENDIX 
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T. R. BRUSSAT 

LOCKHEED -CALIFORNIA CO. 

BEAM THEORY 

+ 

II 

F. ERDOGAN 
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LEHIGH UNIVERSITY 

REISSNER PLATE THEORY 

A 

IV 

GEOMETRIC 

LINEAR 

G. P. ANDERSON/ 
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K. L. DEVRIES 

THIOKOL CORPORATION/ 
UNIVERSITY OF UTAH 

F. E. TASS 

CB 

I 

R. A. EVERETT 
J. D. WHITCOMB 

NASA LANGLEY 

F. E. GAMNAS 

□ 

V 

W. L. HUFFERD 

UNITED TECHNOLOGIES-CSD 

F. E. TEXGAP 

■ 

VI 

C. LOF 

NATIONAL AEROSPACE LAB. 
NLR (THE NETHERLANDS) 


o 
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NONLINEAR 
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INDIAN INSTITUTE OF 
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F. E. GAMNAS (EAST) 

O 
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F. E. VISTA 

fS 

VI 

G. E. LAW 

GENERAL DYNAMICS/FW 

F. E. NASTRAN 

• 
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R. A. EVERETT 
J. D. WHITCOMB 

NASA LANGLEY 

F. E. GAMNAS 

O 

V 

MIXED 
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U. OF MISSOURI-ROLLA 

BEAM THEORY + LINEAR F. E. 
GEOMETRIC NONLINEAR 

a 

IX 
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• Adhesive bondline is 0.13mm thick 


Fiaure 1 - Eaual thickness adherend snerimen 



• Strap adherend is 6.35mm thick v(0,0,z)=0 

•Adhesive bondline is 0.13mm thick 


Figure 2 - Unequal thickness adherend specimen 



EQUAL THICKNESS ADHERENDS 



DEBOND LENGTH, mm 

Figure 3 - Predicted total strain energy release rate versus debond length for 
equal thickness adherend specimens. 



EQUAL THICKNESS ADHERENDS 



DEBOND LENGTH, mm 

Figure 4 - Predicted ratio of mode I to mode II versus debond length for equal 
thickness adherend specimens. 
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Figure 5 - Predicted total strain energy release rate versus debond length for 
unequal thickness adherend specimens. 




UNEQUAL THICKNESS ADHERENDS 



DEBOND LENGTH, mm 

Figure 6 - Predicted ratio of mode I to mode II versus debond length for 
unequal thickness adherend specimens. 




APPENDIX I 


G. P. Anderson and L. P. Abrahamson 
Thiokol Corporation 
Brigham City, Utah 

and 

K. L. DeVries 

Department of Mechanical and Industrial Engineering 
University of Utah 
Salt Lake City, Utah 


APPROACH 

The analysis of the equal thickness adherend cracked-lap-shear specimen 
(CLS-A) was completed using a linear elastic finite element computer program 
(TASS - generated by Morton Thiokol, Inc.)* No special crack tip elements 
were used. A thin row of elements was input at the center of the adhesive; a 
crack was simulated by replacing these elements with "void" elements (that is, 
a zero or near zero stiffness element). 

Mode I and mode II energy release rates (G j and Gjj) were calculated 
using the modified crack closure method outlined in the reference below. The 
grid network consisted of 1,710 quadrilateral elements, three elements through 
the thickness of each adherend and seven through the adhesive. Each quadri- 
lateral element was automatically divided into four linear displacement ele- 
ments by the computer to calculate grid displacements. The grid network as 
shown in Fig. 1-1 contained 0.003 in. x 0.0008 in. quadrilateral elements near 
the crack tip. 

The modified crack closure method required the two equations: 
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where F x and Fy are the forces required to close the crack a distance 
Aa and u x -j, Uyj are the x and y crack opening displacements a dis- 
tance Aa behind the crack tip [1-1]. 

RESULTS 

The resulting energy release rates for five crack depths are presented in 
Table 1-1. 

The initial analysis for the 101.6 mm (4 in.) crack depth used a linear 
displacement element and the grid shown in Fig. 1-2. A total energy release 
rate of 375 J/m^ was obtained. It was later determined that the grid network 
between x = 50.8 and x = 139.7 mm was too coarse to provide an adequate 
beam bending analysis. The proper energy release rate (213 J/m^) was obtained 
by using a quadratic displacement element with a coarse grid similar to 
Fig. 1-2 or by using a finer grid (Fig. 1-1) with the linear displacement 
element. The grid in Fig. 1-2 used 1,710 elements while the grid in Fig. 1-1 
required 2,347 elements. The quadratic displacement element grid used 640 
eight-node elements. 
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TABLE 1-1 


Calculated Strain Energy Release Rates 
With Debond In Middle Of Adhesive 

J/m^ 


Debond 

Length 

in. 

mm* 

0.10 

2.54 

0.25 

6.35 

1.0 

25.40 

2.0 

50.80 

4.0 

101.60 

CLSA 






G I 

35 

42 

57 

95 

48 

G II 

144 

153 

180 

234 

165 

g t 

179 

195 

237 

329 

213 
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APPENDIX II 


T. R. Brussat 

Lockheed-California Company 
Burbank, CA 91520 

APPROACH 

Using elementary beam theory closed-form analytical solutions were 
obtained for the mode I and mode II components of the strain-energy- release 
rate for both CLS specimen geometries. 

Reference II-l, which first introduced the CLS specimen, provides most of 
the equations used here. The assumption is made in Ref. II-l that the length 
of the specimen and crack are large compared to the thicknesses; consequently 
the equations are all independent of crack length. 

The adhesive layer is relatively flexible and relatively thin. The con- 
tribution of adhesive stiffness can therefore be neglected. However, the 
bondline thickness does significantly affect the offset distance characteriz- 
ing the crack-tip eccentricity, and it also significantly affects the moment 
of inertia on the uncracked end of the specimen. Therefore, bondline thick- 
ness is considered in calculating these quantities. 

Analytical Equations 

An exact expression for total strain-energy-release rate of an infinite- 
length CLS specimen is derived in Ref. II-l. The resulting Equation is 

G ° 2 , 

“T ‘ 2b N («)2 

where P is the applied load; b N is the specimen width (measured at the 
bondline); (EA )2 is the tensile rigidity of the strap; and (EA) 0 is the total 
tensile rigidity (lap + strap). 


1 - 


(EA), 

teat; 


(II-l) 
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Ordinary beam theory is used in Ref. II-l to obtain expressions for the 
mode I crack opening displacements and the limiting value of the internal 
bending moment at the crack tip for the infinite-length specimen. The deri- 
vation presented in Ref. II-l could have been modified to satisfy the round- 
robin boundary conditions (i.e. specimen length and loading conditions), but 
for simplicity this was not attempted. The expressions given in Ref. II-l are 
given in terms of two dimensionless parameters, Y 2 and y 0 , which are 
related as follows to the bending rigidities (El ) 2 and (El ) 0 of the strap 
section and the combined (lap + strap) section, respectively: 

x 2 = /p/(EI ) 2 

X 0 * /p/(EI) 0 (1 1—2) 

The mode I crack opening displacement a distance x from the crack tip is 
given by 


y' 


*2 


1 + ( X 




= e 


-A 2 x 



(II-3) 


where y 2 is the centroid location of the strap section, and y Q is the 
centroid location of the combined section. The limiting value of internal 
bending moment at the crack tip is 


(y 2 - y 0 )p 

M = — — 

° 1 + 


(1 1—4) 


An approximation for the mode I component of strain-energy-release rate 
is derived in Ref. II-l under the following assumptions: 
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(1) Gj for a CLS specimen in tension is equal to Gj for a CLS specimen 
subjected to end moments of magnitude M 0 . 

(2) Gjj/Gj = 4/7 for the CLS specimen in pure bending. (This is the 
exact beam-theory result for equal thickness adherends and a zero-thickness 
adhesive layer). 

In accordance with these two assumptions, the Equation for Gj given in 
Ref. II-l is 


G 


2M o 

I = 7b N (EI) 2 


(ei) 2 


(II-5) 


Recently the author has re-examined the second assumption above. For 
purposes of the round-robin analysis study, a second way of estimating Gj was 
proposed, based on the following alternative assumption: 

The problem of the CLS specimen in pure bending can be separated into the 
approximately Pure mode I and approximately Pure mode II problems shown in 
Fig. II-l. 

The pure bending of the CLS specimen under end moments, M 0 , is shown in 
Fig. II-l(a). In (b) the moment M 0 is balanced by a pair of moments, M^ 
applied to the lap and M 2 applied to the strap. The magnitudes of 
and M 2 are each proportional to the ratio of the bending rigidity divided by 
the centroidal distance from the bondline, (El )^/ |y^ | and (EI^/^I respec- 
tively. This creates equal and opposite bending strains along the two faces 
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of the crack, which is an antisymmetric condition. The magnitudes of Mj 
and Mg are given by 


, V^VIyjl 

1 " WTJW[\ + (ED 2 /|y7T 

(II-6) 

M 0 (EDi/|y 2 | 

2 " TETi' 1 /|y 1 | * (ED 2 /|y 2 t 


In Fig. II-l(c), an approximately pure mode I case is shown such that 
superposition of (b) and (c) would lead to (a). Thus, the solution Gj for the 
case shown in Fig. II-l(c) is an alternative approximate solution for Gj for 
the CLS specimen loaded in tension. Tada, et al., (Ref. 1 1-2 ) give the 
following solution for the configuration shown in Fig. II-l(c): 


G 


I 


2b, 



( I 1—7 ) 


In the results that follow. Method 1 uses Eq. (1 1-5) to estimate Gj, 
while Method 2 uses Eqs. (1 1-6) and (1 1-7 ) . 

Gj i j vanishes in this 2-dimensional analysis. Therefore, the mode II 
strain-energy-release rate can be estimated by subtraction: 

6 II = g T - G j (1 1-8) 


The relative sliding-mode displacements, ax, between the crack surfaces 
result additively from the tensile load and bending moment in the strap; the 
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lap is assumed to be stress-free. From the beam theory analysis results given 
in Ref. II-l it can easily be shown that 

Pv - AoX 

^ = + 2(El ) 2 a 2 (1 " e ) (II " 9) 

where h£ is the thickness of the strap. 

RESULTS AND DISCUSSION 

Table II-l gives the strain-energy-release rate components for CLS A and 
CLS B. 

Note that the two different methods of estimating Gj give very different 
results for the unequal adherend case, CLS B. Eq. ( 1 1 -5 ) (Method 1) results 
in a Gj/Gji ratio of 0.36, whereas Eq. (I I -7 ) ( Method 2) gives Gj /G j j = 0.54. 

Since the supporting assumptions for Method 2 seem more valid, but the 
results for Method 1 seem more likely to be correct, it is of interest to sub- 
mit both for the round-robin study. 
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TABLE II-l 


Calculated Values of Strain Energy Release Rate 

J/m^ 



Method I 

CLSA 

Method II 

Method I 

CLSB 

Method II 

G I 

44 

44 

74 


97 

G II 

164 

164 

204 


181 

g t 

208 

208 

278 


278 

G I /G II 

0.27 

0.27 

0.36 


0.54 
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(a) CLS in Pure Bending 




(c) Approximately Pure Mode I 


Figure 1 1 -1 . Approximately symmetric and antisymmetric components for the 
CLS specimen bending. 
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APPENDIX III 


B. Dattaguru and P. D. Mangalgiri 
Department of Aerospace Engineering 
Indian Institute of Science 
Banagalore, India 

APPROACH 

A Finite Element Analysis was carried out on a two dimensional plane 
strain idealization of the cross-section of these joints. The eccentricity of 
load transfer in these joints causes large rotations and so a geometrically 
nonlinear finite element analysis was employed (Refs. III-l, III-2, and 
I I 1-3). The basic approach used was Lagrangian where the displacements are 
referred to the underformed configuration of the structure. The geometric 
nonlinearity was introduced in the strain-displacement relations as 



y - l u _ + 9v + 9u 3u 3y 9v 

'xy 9y 9x 3x ay ax ay 

where u and v are displacements along x and y axes. Isoparametric 
elements were used. The tangent stiffness matrix at any stage of deformation 
was calculated as 


k t = 


K o + 


k l + 


(III-2) 
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where K Q , K|_ and Kg are the linear, large displacement and geometric 
stiffness matrices. A Newton-Raphson iterative scheme was employed and the 
tangent stiffness metrix was updated after every four iterations. 

CALCULATIONS OF Gj AND G n 

The mode I and mode II strain-energy-release rates Gj and Gjj were calcu- 
lated based on virtual crack extension method (Ref. 1 1 1-4). In order to 
account for large rotation of the debond, the components of forces and dis- 
placements in the directions along and normal to the center line of the 
deformed debond configuration were used to calculate Gj and Gjj (Fig. III-l). 
Thus 


8 .Ip ^ - V 3 l 

I 2 y' b a a 


'II 


1 P 
7 P x' 


(uj - U3 ) 
b a a 


( 1 1 1 -3 ) 


where u' and v 1 are displacements along x 1 and y' axes as shown in 
Fig. III-l. 


Finite Element Model 

A typical finite element mesh used for the present analysis is shown in 
Fig. 1 1 1-2. The mesh has 371 nodes and 320 elements. Other models varied 
between 305 to 375 nodes and 258 to 320 elements. In all the models, the thin 
adhesive was divided into two layers of elements across the thickness. The 
applied loading was assumed to be uniform on the loaded end and was distri- 
buted as a consistent load vector at the nodes. 
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NUMERICAL RESULTS 


Analysis was carried out for debond lengths a = 0.1, 0.25, 1.0 and 
4.0 in. The strain-energy-release rates for all these cases are shown in 
Table III-l. 
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TABLE III-l 


Calculated Strain Energy Release Rates 
With Debond In Middle Of Adhesive 

J/m^ 


Debond 

Length, 

in. 

mm. 

0.10 

2.54 

0.25 

6.35 

1.0 

25.40 

4.0 

101.60 

CLSA 





G I 

37 

37 

37 

38 

G II 

150 

150 

151 

153 

g t 

187 

187 

188 

191 


CLSB 


G I 

49 

49 

50 

44 

G II 

224 

225 

227 

218 

g t 

273 

274 

277 

262 


37 



CO 

00 




Figure III-l. Transformed coordinate system for 6 calculations 




APPENDIX IV 


F. Erdogen and P. Joseph 
Lehigh University, Bethlehem, PA 18015 

ASSUMPTIONS 

The stress analysis of the round robin cracked-lap-shear specimens is 
solved under the following assumptions: 

(a) The adherends are approximated by Reissner plates. That is a plate 
theory taking into account the transverse shear effects rather than continuum 
elasticity is used in formulating the problem. 

(b) The problem is assumed to be one of plane strain; that is e z is 
assumed to be zero for the entire specimen. 

(c) The adhesive is assumed to be an elastic layer in which the thickness 
variation of stresses is neglected. In formulating the adhesive a slight im- 
provement is made over the conventional uncoupled tension-shear spring model 
by taking into account the effect of the average in-plane strain e x . 

Partial reasons for adopting the particular analytical model for the 
adhesive joint are as follows: 

(1) Generally, in adhesively bonded structures, the thickness of the ad- 
hesive is very small compared to the thicknesses of the adherends which, in 
turn, are small compared to the in-plane dimensions of the joint. The "plate" 
theory is known to deal quite satisfactorily with structures having such 
geometries. 

(2) The plane problem can be solved in closed form (Ref. IV-1). 

(3) The technique can be extended to treat adhesively bonded joints with 
complicated geometries and to take into consideration such effects as viscoe- 
lastic behavior of the adhesive (Ref. IV-2). 
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(4) Even though the model is not suitable for the calculation of a 
"stress intensity factor", it is suitable for the calculation of the "strain- 
energy-release rate". 

(5) For one specimen geometry tested the tensile and shear stresses in 
the adhesive obtained from the plate model appear to be in good agreement with 
those given by the finite element method (Ref. IV-1). 


APPROACH 

The two adherends are assumed to be "plates" under in-plane deformations 
and bending. The equal ibri urn equations for the lap and the strap adherend may 
be expressed as follows: 




dM lx 

T5T =( >1* 


h l + h o 
2 


T 


(iv-i) 


dN 


2x 


dQ 


dx 


= -t, 


2x 


dM 


dx 


= -a. 


2x 


dx 


= Q 


h 2 + 


2x 


(IV-2) 


where for i = 1 and i = 2 N ix , Q ix , M ix , are the stress and moment resul- 
tants in the lap and strap adherends, respectively, and t and a are the 
shear and the normal stress in the adhesive. The stress and moment resultants 
are related to the x, y-components of the displacements u^, v^ and to the 
rotations B ix , (i = 1,2) by 


du, dv. 

dx Vix* dx 


+ 8. 
P 1X 


G 1x d6 fx 
B, ’ dx 


D lV < W ’ 2) 


(IV-3) 


C i * 


1- V • 

v ix 1Z 


B r5 h t G ixy = 


12 ( 1- v* v • ) 
1 IX iz ; 


^i7 


(IV-4) 
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Assuming that y-dependence of the strains e x , Ey , and y xy 1n the 
adhesive is negligible, from kinematical considerations it may be shown that 

V 1 “ v 2 1 1 
e y FT ’ Y xy " (u l " 7 h l 6 r ' u 2 “ 7 h 2 B 2x^ /h o’ 

^ A 


e 


x 


_ , du l h l dB lx + du 2 , 

■ ' t sr + ir + 


h 2 dB 2x 
2 dx 


)/ 2 , 


( IV— 5 ) 


where h Q , h^ and hg are the thicknesses of the adhesive, the lap adherend 
and the strap adherend, respectively. If E, and v denote the elastic 
constants of the adhesive, its stress-strain relations may be expressed as 


2 


l-v-2v 



e , Y 
x’ 'xy 


2(1+V) 
E 


(IV-6) 


By simple eliminations. Equations (I V-l ) - ( IV-6 ) may be reduced to a 
system of differential equations for the functions x(x) and o(x) which can 
then be solved in closed form (see Ref. I V-l for details). 

Strain Energy Release Rate 

In an elastic structure containing a flaw of "area" A, ignoring the 
dynamic effects, the energy balance Equation may be expressed as 

G = ^ (U-V) = r F 

where U is the work done by the external forces, V is the stored elastic 
energy and yf 1S the fracture energy of the material. In Equation (IV-7) 
the left-hand side represents the externally added or internally released 
energy available for fracture, and yp is the measure of the fracture resis- 
tance of the material. If the bulk of the structure undergoes elastic 
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deformations, it is known that G is the same under "fixed grip" and "fixed 
load" conditions. Thus, G can be calculated similar to the crack closure 
energy by considering the advance of the debond front, and by assuming fixed 
grip conditions. 

As the debond front advances by a length da, dl) = 0 and dV (per unit 
crack front) may be calculated by relaxing the stress state in the adhesive 
for a volume h 0 (1.0) (da) and the surface tractions a(X) and x(X) acting 
on the adherends along the debond area da to zero. The strain energy released 
by the adherends due to the relaxation of the tractions a and t may be 
expressed as 

dVj = - £ ada (« tl + 6 t2 ) xda (« $1 + 6 g2 ) (IV-8) 

where 6^.^ and 6 S1 -, (i=l,2) are the y and x-components of the displacements 
of the adherend surfaces at the debond region da due to the removal of the 
tractions ada and xda and the minus sign is due to the fact that during 
the release process, the directions of the forces and the displacements are 
opposite to each other. 

The strain energy released from the relaxation to the adhesive may be 
obtained from 

a+da h rt 1 

dV 2 = - / dx / 0 dy / dz W (I V-9 ) 

a oo 

where W is the strain energy density in the adhesive. In the model used and 
for the plane strain problem under consideration, the adhesive stresses are 
assumed to be independent of y and z. Thus, Eq. (IV-9) can be expressed as 

dv 2 = -Wh 0 da = - \ (o xtx ♦ < IV - 10) 
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where W is calculated at the debond front. For the problem under considera- 
tion we have (REF. IV-1): 


dV 


2 



(a) 


1-v 


o^(a) l-v-2v^ 

E T=T + 




( IV— 11 ) 


where E, v and G are the elastic constants of the adhesive and e x , a, 
and t are calculated at the debond front a. 

We now observe that the plate theory would give the displacements in the 
adherends as follows: 


°da , xda (1-v-j ) (1 =1.2) 

6 ti “ E7" 6 si = — r — 

i i 


( IV— 12 ) 


From Eqs. ( I V— 8 ) and (IV-12) it then follows that as da+o the strain-energy- 
release rate dVj/da contributed by the adherends would approach zero. Since 
dV = dV 1 + dV 2 and dA = da (per unit debond front), from Eqs. (IV-7) and (IV- 
11) we obtain 


G = 


da 


(U-V) = 


o 

T~ 


~T7 


l-v-2v^ o' ( a ) 


x (a) 2(1 
E 


=3 


( IV- 13) 


The strain-energy-release rate G calculated from Eq. (IV-13) would be 
equivalent to the conventional Gj + Gjj. It should be noted that if at the 
debond front the adhesive is in compression (which generally is the case if 
the bending stiffness of the lap adherend is greater than that of the strap 
adherend and if there is no transverse constraint in the strap adherend), then 
Eq. (IV-13) should be modified as follows: 
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G = G 


ii s -r 


—j^r 


t (a)2(l+v) 


, (o(a) < 0) 


(IV-14) 


In this case the problem is equivalent to Kj = 0 in a crack under mixed mode 
conditions, and the effect of possible crack surface friction is ignored. 


RESULTS 

As pointed out in the previous section, the plate model used in this 
study can give only the total strain-energy-release rate G rather than Gj 
and Gjj separately. Furthermore, since plane strain conditions are assumed to 
prevail in z-di recti on, in the present solution Gjjj = 0. For the trans- 
versely constrained loading condition and dimensions shown in Figs. I V— 1 and 
IV-2 the results for the complete range of the debond length are given in 
Table IV-1 (see also Table IV-2). Here the adherends are assumed to be alumi- 
num (E = 72.450 MPa, v = 0.33) and the elastic constants of the adhesive are 
E = 1932 MPa, v = -.40. To indicate the overall trend, the results are also 
shown in Fig. IV-3. 

From the expression for the strain-energy-release rate given by 
Eq. (IV-13) and from the calculated results, it was observed that the contri- 
bution of the first term (involving e x ) to G is approximately two orders of 
magnitude smaller than that of the remaining terms. Thus, if the effect of 
e x is neglected, the second and third terms in Eq. (IV-13) may be interpreted 
as Gj and G j j , respectively. Partial results giving the individual contribu- 
tions of the terms involving o, t, and e x in Eq. (IV-13) are given in 
Table IV-2 and are labeled as Gj, Gjj, and G e , respectively. Note that G 
is the sum of these three terms. 

Fig. IV-1 also shows the strain-energy-release rate for the cracked lap- 
shear specimen without the transverse end constraint. In this problem the 
transverse shear force Q at the end is zero and the specimen is free to 
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"bend". Consequently, the normal stress in the adhesive is zero for the spec- 
imen with equal thickness adherends and compressive for the specimen with 
unequal thickness adherends. Thus, for these two specimens G is calculated 
from Eq. (IV-14). A peculiar result observed in these calculations was that 
in varying the debond length a from zero to 229 mm, G turned out to be con- 
stant, namely 

G = 45 J/m 2 for CLS.Afhj = h 2 ), 

G = 74 J/m 2 for CLS.B (h 1 = 2h 2 ) 

In the plate model adopted in this study, it is assumed that in the 

"debonded" part of the joint the adhesive layer is completely unloaded. 
Therefore, in this model, the results are not sensitive to the location of the 
"crack" in the adhesive and the calculated "crack opening displacement" does 
not have the conventional meaning. 
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TABLE IV-1 


Calculated Strain Energy Release Rates 
With Debond in Middle of Adhesive 

J/m 2 


Debond 

Length 

in. 

mm. 

- 

0.1 

2.54 

0.25 

6.35 

1.0 

25.40 

4.0 

101.60 

CLSA 

G I 


39 

43 

62 

51 

G II 


145 

151 

180 

163 

G e 


0.3 

0.3 

0.4 

0.4 

g t 


184 

195 

243 

214 

CLSB 

G I 


102 

110 

126 

19 

G II 


284 

294 

311 

175 

G e 


0.6 

0.6 

0.6 

0.3 

g t 


387 

404 

437 

194 
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APPENDIX V 


R. A. Everett, Jr.* and J. D. Whitcomb 
NASA Langley Research Center 
Hampton, VA 23665 

APPROACH 

To analyze the cracked-lap-shear (CLS) bonded joint configurations speci- 
fied in this round robin, a two-dimensional finite element program called 
GAMNAS (Geometric and Material Nonlinear Analysis of Structures) was used. 

This program was developed at NASA Langley to support fracture mechanics 
studies of debonding and delamination and is documented in Refs. V-l and V-2. 
This study used a nonlinear geometric analysis assuming plane strain 
conditions. To calculate strain-energy-release rates, GAMNAS uses a crack 
closure technique like that reported in Ref. V- 3. This is done by using the 
forces transmitted through the node at the crack tip and the relative dis- 
placements of the two nodes on the crack boundary closest to the crack tip to 
calculate the energy required to close the crack. 

No special crack tip elements are used in GAMNAS. For this analysis a 
4-node isoparametric quadrilateral element was used. This finite element 
program has options for full and selective reduced integration. In this anal- 
ysis selective reduced integration was used to improve the element's perfor- 
mance in modeling bending type deformations. 

The mesh for the CLSB specimen with a 101.6 mm debond is shown in 
Fig. V-l. All the analyzed configurations were modeled similarly. All models 
had about 1000 elements with the thickest adherend having 9 elements through- 
the thickness and the thinnest having 7. The adhesive had 4 elements through- 
the thickness. At the debond tip the elements had an aspect ratio of one with 
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the length of the element being 0.032 mm. The debond was modeled as a 
perfectly smooth crack between the second and third elements in the adhesive 
layer (in the middle of the adhesive). 

RESULTS 

The strain-energy-release rates calculated in this analysis are presented 
in Table V-l. Both linear and nonlinear results are given for mode I, mode 
II, and the total strain-energy-release rates at the four debond lengths 
analysed. 

The most significant observation from the results in Table V-l is that 
the Gy calculated from the nonlinear geometric analysis is almost constant 
with debond length, whereas, the linear results show Gy to vary with debond 
length with a maximum value at one inch. The ratio of G j /G j j behaves in a 
similar manner. In general, for both configurations the nonlinear value of G 
is less than the linear value. The results also show that the CLSB configura- 
tion with the thicker lap adherend gives a higher value of Gy. 
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TABLE V-l 


Calculated Strain-Energy-Release Rate With Debond 
In Center Of Adhesive, J/nr. Both Geometic 
Linear And Nonlinear Results Are Presented. 


Debond 


length 

in. 

mm. 

0 . 

2. 

linear 

.10 

54 

nonlin. 

0 . 

6 

1 i near 

.25 

35 

nonlin. 

1 . 

25 

linear 

.0 

4 

nonlin. 

4, 

101 

linear 

.0 

6 

nonl i n . 

CLSA 









G I 

37 

39 

44 

39 

60 

39 

47 

40 

G II 

145 

147 

149 

149 

180 

151 

163 

152 

g t 

182 

186 

193 

187 

240 

189 

212 

193 

CLSB 









G I 

100 

51 

110 

51 

123 

53 

19 

46 

G II 

284 

221 

291 

222 

312 

226 

175 

217 

g t 

383 

271 

401 

273 

434 

279 

194 

263 


51 
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Figure V-l. 


Finite element mesh for 101.6 ntn debond in the unsymmetri 
adherends specimen configuration. 




APPENDIX VI 


William L. Hufferd 
United Technologies 
Chemical Systems Division 
San Jose, California 

APPROACH 

Linear elastic analyses of the two cracked lap shear geometries were 
conducted using two different finite element codes: TEXGAP (Ref. VI-1) and 

VISTA (Ref. VI-2). Geometrically nonlinear analyses were conducted for one 
cracked lap shear geometry using VISTA. 

The version of TEXGAP used at CSD calculates stress intensity factors in 
one of three ways: (1) using a hybrid crack element, (2) using contour inte- 

gration, and (3) using a singular crack element. The hybrid crack element was 
used in the current calculations. This element is based on a displacement 
formulation in which the displacements are interpolated over the boundary of 
the element and the stresses are interpolated over the interior of the ele- 
ment. Mode I and mode II stress intensity factors are calculated directly. 

The hybrid element is an 11-node, square element with its local coordi- 
nate system located at the center of the element (i.e., the crack tip). The 
element may be used for plane stress or plane strain geometries or axisym- 
metric geometries at large radius. It has been reformulated for 
incompressible or nearly incompressible materials, and it includes thermal 
loadings, but excludes body forces. The crack surfaces within the element are 
assumed traction free. 

The hybrid crack element used 15 interpolating, quadratic displacement 
functions for boundary displacements which insure interelement compatability 
with adjacent quadratic elements, and 19 different interpolating stress func- 
tions for the interior of the element which identically satisfy both 
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compatability and equilibrium. Symmetric and antisymmetric stress distribu- 
tions are included as well as mode I and mode II stresses. The stress fields 
model the square-root singularity and also incorporate eigenvalues greater 
than one. The angular dependence of the stress distribution is also appro- 
priately modeled. 

VISTA is a finite element code for the solution of two dimensional (axi- 
symmetric or generalized plane strain) quasistatic viscoelastic stress anal- 
ysis problems with small strains and small or large displacements. The 
singular element used for fracture mechanics analysis with VISTA is based on 
the interpolating shape functions given by Stern (Ref. VI-3). It is a six- 
node subparametric triangle. The corner nodes are used for geometry inter- 
polation and the midside and corner nodes are used for displacement inter- 
polation. Thus, it is a straight sided element. The order of the singularity 
(one-half for linear elastic, isotropic materials) is input by the user. The 
displacements are quadratic along the side of the element opposite the singu- 
lar point, thus conforming with the quadratic isoparametric elements in VISTA. 

The implementation of the element in VISTA uses the standard element 
stiffness and load vector routines. Thus, the constitutive properties are 
handled by the same routines as regular elements so that any material can be 
used in the singular elements. The numerical integrations in the singular 
element use a specially derived quadrature rule. 

The code outputs the coefficients of the displacement field; i.e.. 



where X is the order of the (user specified) singularity and p and 0 are 
polar triangular coordinate. For the mixed mode problem in an isotropic 
linearly elastic solid with X = 1/2, the plane strain near field 
displacements are given by: 


K = irjfa cos l[ x " 2v + sin2 Z_ 


+ IT* 


Irjh sin l[ : 


2 - 2v + cos 


?] 


U y = Trjfc C0S I [ 2 ' 2v + s1 " 2 l] 

* -¥-Jh cos 1 [- 1 + 2v + si " 2 l] 


where G is the shear modulus and v is Poisson's ratio for the material. 
Selecting the crack faces at 0 = ±ir for evaluating Kj and Kjj leads to 
the simple expressions: 

K. = K-fa f (0) 

1 4(1 -v 2 )V v 

k„ = — u f (0) 

11 4(l-v 2 yV u 
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These latter expressions were used to calculate stress intensity factors from 
VISTA output for f u (0) and f v (0)» from which mode I and mode II strain- 
energy-release rates, Gj and Gjj, respectively were calculated using: 

(1 - V 2 ) K 2 

G I E 1 

c . d-v 2 )^ 

G II E 

DESCRIPTION OF ANALYSIS CONDUCTED 

A typical deformed finite element model is shown in Fig. VI-1 for the 
CLS-B geometry with a 25.40 mm cohesive crack through the midplane of the 
adhesive layer. The model contains 397 elements. Four elements were used 

through the 5-mil thickness of the adhesive layer. The mesh in the neighbor- 
hood of the crack tip is shown in Fig. VI-2. The crack-tip region itself was 
modeled with eight singular triangular elements. The singular elements were 
overlaid with the hybrid crack element as shown in Fig. VI-3 for the TEXGAP 
analyses. All other finite element models for analyses of both crack lap 
shear geometries: CLS-A(3.18 mm thick lap and strap) and CLS-B (6.35 mm thick 

lap and 3.18 mm thick strap) were similar, except that the 4x4 fine mesh was 
moved with the crack tip for other crack lengths. 

Linear elastic analyses were conducted for both CLS geometries for 2.54, 
6.35, 25.4 and 101. 6 mm cohesive cracks in the adhesive. A geometrically 
nonlinear analysis was conducted for a 25.4 mm crack in the 6.35 mm lap and 
3.18 mm strap CLS geometry (CLS-B). All analyses assumed plane strain 
conditions. 
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DISCUSSION OF RESULTS 


One linear elastic analysis was conducted using VISTA to provide a 
baseline for comparison with TEXGAP analysis results. The two codes gave 
virtually identical stresses and displacements and the stress intensity 
factors calculated from the coefficients of the displacements from VISTA were 
within 2 percent of those obtained from the hybrid crack element used with the 
TEXGAP analyses. 

Table VI-1 summarizes the calculated strain-energy-release rates. A 
slight maximum is observed at a crack length of about 25.4 mm. The thicker 
lap of CLS-B results in higher mode I and mode II values than is observed for 
the equal thickness adherends. 

A geometrically nonlinear elastic analysis was conducted using VISTA in 
which the total load was applied in two load steps. Convergence for the first 
load step took six iterations, while that for the second took four iterations. 
A major effect of the nonlinear analysis was to "smooth" the distortion of the 
singular crack-tip elements on the free surface of the crack face. As a 
result, Kj dropped approximately 50 percent from the VISTA linear analysis 
while Kj j changed only about 20 percent. 

As a final remark, the absolute values of the stress intensities computed 
from the nonlinear analysis should be carefully interpreted. These calcula- 
tions were made assuming that the order of the crack-tip singularity was one- 
half, the same as for a linear elastic analysis; and the calculations were 
made using the same displacement equations as used in a linear elastic 
analysis. The validity of these assumptions for this nonlinear problem is not 
known. In general, in this situation, more reliable results would be obtained 
from a patch-independent integral calculation, such as a J-integral, which has 
demonstrated validity for nonlinear elasticity problems. 
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TABLE VI-1 


Results With Debond In Middle Of Adhesive Using TEXGAP 

(Geometric Linear) 

J/m 2 


Debond 

Length, 


in. 

mm. 

0.10 

2.54 

0.25 

6.35 

1.0 

25.40 


4.0 

101.60 

CLSA 

G I 

35 

37 

60 


42 

G II 

142 

149 

172 


156 

g t 

177 

186 

232 


198 

CLSB 

G I 

88 

114 

116 

47* 

12 

G II 

282 

298 

310 

207* 

168 

g t 

370 

412 

426 

254* 

180 

^calculated using VISTA (Geometic 

Nonlinear) 

at a=25.40 mm. 
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Figure VI 1. Deformed Finite Element Model of CLS-B (6.35 mm Lap and 3.18 mm Strap) 
Geometry with 25.4 mm Cohesive Crack in Adhesive Layer (Deformed Scale 
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Figure VI-2. Crack-Tip Finite Element Model 





Figure VI-3. Details of the hybrid crack element overlaid singular 
triangular elements. 


APPENDIX VII 


George E. Law 
General Dynamics 
Fort Worth Division 
Fort Worth, Texas 

APPROACH 

Non-linear geometric finite element analysis was performed with MSC/ 
NASTRAN using "Solution 64". This is a non-linear geometry solution sequence 
for large displacement/small strain applications. The solution technique is 
based on the geometric stiffness approach. Details of this method are con- 
tained in Section 2.9 of the "MSC/NASTRAN Application Manual Volume 1", The 
MacNeal-Schwendler Corporation, Los Angeles, California, May 1983. 

The two-dimensional finite-element model constructed for the analysis of 
the cracked lap shear specimen analysis is shown in Fig. VII-1. Isoparametric 
3-node triangles and 4-node quadrilateral elements were employed in this 
model. The model was constructed using the load and out-of-plane directions 
to describe the 2D space. A state of plane stress was assumed. Fig. VII-1 
shows the outline of the exterior surface of the model accentuating the crack 
line. A blow-up of the mesh at the crack tip area is also shown. Two models 
were made: one with a short, incomplete glue line and the second in which the 
complete glue line was modeled downstream from the crack tip. The elements in 
the crack-tip region were 0.0159 x 0.0159 mm. 

The boundary conditions applied to the model were 
u(0,y) = 0 
v(0,0) = 0 
v (L ,0 ) = 0 

The loading was applied as a tensile force of 2268 kg. (2500 lb.) at x = L; 
the force was uniformly distributed through the thickness of the model. 
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The modified crack closure technique (Ref. VII-1) was used to calculate 
the mode I and mode II components at the crack tip. Referring to Fig. VI 1-2 , 
the nodal forces at "f" and the displacements at the first upstream nodes, "g" 
and "h" , are combined to calculate the work to close the crack. The mode I 
and mode II components of the energy release rate are calculated as 

Gj = Fy dv/2Aa 
Gji = F x du/2Aa 

where 

dv = v„ - v/ 

9 h 

du = u - u/ 
g h 

F x and F y are the forces in the respective X and Y directions that 
resist the crack against opening, and Aa is the distance between node "f" 
and nodes "g" and "h" in the undeformed state. 

In the case of the non-linear geometric analysis, the X- and Y- axes for 
the mode I and mode II components must be defined in the deformed state. 

Fig. VI -3 shows the relations required to determine a new X' -Y' coordinate 
system aligned with the crack. The results of the finite element analysis are 
transformed into the X' -Y' coordinate system through standard tensor trans- 
formations and the mode I and mode II components of the energy release rate 
are then calculated. 

RESULTS 

Two glue line models were used in this analysis. The short glue line 
model was found to produce erroneous results when compared with the full glue 
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line model. This short glue line model was exercised for all of the analysis 
cases of the round robin. Those results, although erroneous in magnitude, are 
presented in Table VII-1 to describe trends in the CLS analysis. The full 
glue line model was only exercised for two cases of the round robin. Since 
the results for that analysis agreed with results of other round-robin partic- 
ipants using nonlinear geometric analysis, it is inferred that all of the 
other cases would also agree. 

Table I presents the energy release rates calculated for the equal thick- 
ness and unequal thickness adherend configurations. These values were calcu- 
lated from the short glue line model and are only of value in that the trend 
of the analysis is represented. Two observations can be made from these 
results: the unequal thickness adherend shows greater mode II than the equal 

thickness adherend coupon, and the calculated energy release rates are essen- 
tially constant with crack length except for a slight perturbation in the 
short crack length range. The increase in the mode II component with increase 
in the thickness of the lap of the coupon is a direct result of increased 
stiffness (EA) of the lap. The increase in lap stiffness causes greater shear 
transfer across the bondline. 

The relationship between energy release rate and crack length can be de- 
scribed based on concepts of self-similar crack propagation. For long crack 
lengths, the highly stressed zone-surrounding the crack tip does not interact 
with the end boundaries. Also, since the crack is parallel to the load di- 
rection, the net section is not reduced. Thus, it is inferred that the energy 
release rate should be independent of crack length. For the short crack 
lengths considered, the length of the crack is on the same magnitude as the 
thickness of the adherends. In this case, interactions between the stresses 
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around the crack tip and the boundary conditions, particularly the free 
boundary at the end of the lap, can be anticipated. 

Table VII-2 gives the full glue line model results obtained for the con- 
stant thickness adherend case in the short (2.54 mm) and long (101.6 mm) crack 
length configurations. For the short and long crack length cases, Figs. VI 1-4 
and VI-5 respectively show the relationship between the calculated energy 
release rates and the number of iterations in the nonlinear solution. Itera- 
tion number 1 represents the linear solution and numbers 2, 3, 4, and 5 repre- 
sent each nonlinear iteration. These figures show that the solution for the 
energy release rate converges after two nonlinear iterations demonstrating 
that a linear solution is invalid. 

Conclusions of this study are that the energy release rate is relatively 
constant with crack length in the cracked-lap-shear specimen. Also, the full 
glue line must be modeled in the adhesive CLS coupon to obtain proper internal 
shear transfer in the analysis. 
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TABLE VII-1 


Calculated Strain Energy Release Rates For Short Bondline* 
With Debond In Middle Of Adhesive 
(*Use for analysis trends only) 

J/m 2 


Debond 

Length 

in. 

mm. 

0.10 

2.54 

0.25 

6.35 

1.0 

25.40 

4.0 

101.60 

CLSA 

G I 

35 

36 

36 

37 

G II 

83 

84 

85 

6 

g t 

118 

120 

121 

123 

CLSB 

G I 

36 

38 

37 

33 

G II 

124 

128 

127 

124 

g t 

160 

166 

164 

157 
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TABLE VI 1-2 


Calculated Strain Energy Release Rates For Full Glue Line 
With Debond In Middle Of Adhesive 

J/m^ 


Debond 
Length 
i n. 
mm. 

0.10 

2.54 

0.25 

6.35 

1.0 

25.40 

4.0 

101.60 

CLSA 

G I 

45 

_ — 


47 

G II 

176 

— 

— 

181 

g t 

221 

-- 

““ 

228 
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Figure VII-1. Finite element models for the cracked lap shear specimen. 
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Figure VII-2. Mesh at the crack tip. 




Figure VII-3. Nonlinear analysis requi 
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Figure VI 1-4. Convergence study for full glue line r lel-crack length 
a = 2.54 mm. 
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Figure VII-5. Convergence study for full glue line model-crack length 
a = 101.6 mm. 
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APPROACH 

A finite element program, ASKA (Ref. VIII-1) using three-dimensional 
strain elements (isoparametric) was used. K- and G- values are assumed to be 
related to displacements in the crack tip fields according to linear elastic 
fracture mechanics theory as described by Paris and Si h [VIII-2]. 

At the debond frontline, elements are applied with singular strain field 
at the tips (subnet 1). Other subnets (2 ... 5) consist of normal hexagonal 
27-node elements. The debond size is varied by changing x-coordi nates of the 
subnets in front of or behind the debond tip. The mesh is shown in 

Fig. VIII-1. Gj, Gjj, and Gjjj are derived from K-values by Eq. (VIII-1) for 
mode I 


G 


I 




(VIII-1) 


K 1 is derived from displacements of nodes close to the debond tip, using 
Eq. (VIII-2). 


K r (r) 


(U 1 " U 2 } E^ 

AT 8(1-v 2 ) 


(VIII-2) 


where u lt u 2 are displacements at a distance r from the debond tip at the 
lap- or strap-side, respectively. The Kj-value at the tip is found by linear 
extrapolations of Kj(r)-values, especially from nodes of the crack tip ele- 
ments. Similar formulas [VIII-l, VIII-2] are valid for mode II and mode III 
within the latter case: (1+v) instead of (1-v 2 ). 
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RESULTS 


Table VIII-1: G-values are derived for 3 positions along the debond line 

across the specimens width, (z = 0 is the mid-plane of the bar). The Gj 0 tal" 
value, however, is derived using the virtual crack extension method [VI 1 1-2] 
and the equation: 

G = (VII 1-3 ) 

where dU is the elastic energy variation by a very small local cracksize 
variation dA. 


dU = u [S - S*]u 


where u is the displacement vector and [S - S*] is the variation of the 
stiffness by virtual debond variation. 

This total energy release value is expected to be considerably more accu- 
rate than the separate Gj, Gjj, Gjjj-values obtained. However, the distri- 
bution over three modes can not be found in this way. 

With respect to the total strain-energy-release rate results 
(Table VIII-1) we propose some adaption in order to properly compare these 
three-dimensional results with other two-dimensional data. Therefore, a 
"weighted average" of these results for discrete positions along the debond- 
zone frontline is proposed in the following form: 

G tot(av) = * 4 G tot(z=0) + * 5 G tot(z=.4) + ‘ l G tot(z=.5) 
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Different factors refer to different areas of virtual debond extensions 
(Fig VI I 1-2 ) . The corrected results are given in Table VIII-2. 

A remarkable effect found by the three-dimensional calculation is the 
"closure" of the debond-opening at the edges of the specimen, as seen by 
detail-observations of cross-sectional deformations of both adherends, in a 
plane very close behind the debond tip (Fig. VII 1-3 ) . Local Gj-values 
decrease rapidly from the mid-plane towards the side of the bar, whereas 
Gin-values increase. 

Unfortunately, this phenonema was not studied in more detail, i.e. by 
using finer meshes, or non-linear calculation. 
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TABLE VIII-1 


Calculated Strain Energy Release Rates 
With Debond In Middle Of Adhesive 

J/m^ 


Debond 

Length, 

in. 0.10 0.25 1.0 4.0 

mm. 2.54 6.35 25.40 101.60 


Z/ 

7 1.0 in. 

0 

mi d 

.4 

.5 

side 

0 

mid 

.4 

1 

.5 

side 

0 

mid 

.4 

.5 

side 

0 

mid 

.4 

.5 

side 

CLSA 

G I 

14 

7 

0 

17 

8 

0 

28 

12 

.05 

23 

10 

0 

G II 

119 

105 

122 

135 

113 

136 

154 

131 

155 

141 

115 

136 

G III 

0 

3 

33 

0 

13 

38 

0 

16 

47 

0 

16 

43 

g t 

205 

187 

194 

224 

208 

217 

281 

255 

255 

247 

222 

222 
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TABLE VI 1 1 -2 


Total Strain Energy Release Rate for CLSA 


Debond length 
a, mm 

G tot (averaged) 

J/m^ 

2.54 

1.95 

6.35 

216 

25.40 

266 

101.60 

232 
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Fig. VIII - 3 Cross-sectional deformation 


APPENDIX IX 


S. Mall 

Department of Engineering Mechanics 
University of Missouri 
Roll a, MO 

APPROACH 

The analysis of cracked-lap-shear specimens requires a geometric non- 
linear method [IX-1, IX-2], In the present approach, this nonlinear analysis 
was conducted by combining a simple nonlinear analysis (based on strength-of- 
materials theory) with a linear finite element analysis in the following 
manner [IX-3], Figure IX-l(a) shows the cracked-lap-shear specimen which is 
to be analyzed with a geometric nonlinear method to account for the 
deformation which responds nonlinearly to the applied load geometry. Figure 
IX-1 (b) shows a small region of this specimen near the debond front. This 
region near the debond front can be analyzed with a linear method provided 
moments, axial and shear loads acting on the boundary are obtained from a 
nonlinear analysis. 

A nonlinear analysis based on a simple strength-of-materials theory was 
developed to compute moments, axial and shear loads acting on the boundary of 
the small region of the CLS specimen as shown in Fig. IX-l(b). This nonlinear 
analysis is the extension of a previous analysis [IX-4] . The previous anal- 
ysis was for an infinitely long CLS specimen (i.e. independent of debond 
length), while the present one accounts for the finite length of the specimen 
and the debond. Fig. IX-2 shows the CLS specimen and the deformation of 
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its centroidal axis. From the simple beam theory, the following expressions 
for lateral deflection of centroidal axis were obtained. 

x > 0 


y(x) 


A 2 cothU 2 £ 2 ) + Aq c o t h ( Aq SLq ) 


sinh (Ag^-Agx) 
sinh( a 2 t 2 ) 


x < 0 


y(x) 


-(y 2 -y 0 )A 2 coth(A 2 jj 2 ) 

A 2 C0th(A 2 £ 2 ) + AgCOth ( Ag 


sinh(A 0 i! 0 -A 0 x) 
si nh( Aq£q) 


a 2 -film, 

Aq V p/(EI) o 


Using the above equations, the moments, axial and shear loads acting on the 
small region shown in Fig. IX-l(b) were computed. Thereafter, this region 
with computed boundary loads and moments was analyzed with a two-dimensional 
linear elastic finite element analysis. 

The finite element mesh consisted of 510 four-node, isoparametric quadri- 
lateral elements and had 1200 degrees for freedom. The length of region ana- 
lyzed with FEM was 6.35 mm on each side of debond front (i.e. total length of 
12.7 mm). The analysis was conducted under plane strain condition. The adhe- 
sive was modeled with four layers of elements. The smallest element size near 
the crack tip was 0.0318 x 0.0318 mm. The strai n-energy-release rates Gy, Gj , 
and Gjj in the FEM analysis were computed using a virtual crack closure tech- 
nique (Ref. IV-5). 
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RESULTS 


The computed strain-energy-release rates (Gj, Gj, and G n ) for both CLS A 
and CLS B specimens are presented in Table IX-1. 

DISCUSSION 

A simple approach for analysing the geometric nonlinear problem of the 
CLS specimen is suggested by combining a simple nonlinear analysis (based on 
the strength-of-materials theory) and linear finite element analysis. 
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TABLE IX-1 


Calculated Strain Energy Release Rates With 
Debond In Middle Of Adhesive 

(J/m 2 ) 


Debond 

Length, 

in. 

mm. 

0.10 

2.54 

0.25 

6.35 

1.0 

25.40 

2.0 

50.80 

4.0 

101.60 

CLSA 

G I 

42 

43 

44 

44 

45 

G II 

149 

149 

151 

151 

152 

g t 

191 

192 

195 

195 

197 

CLSB 

G I 

51 

52 

54 

57 

68 

G II 

200 

202 

205 

209 

222 

g t 

251 

254 

259 

266 

290 
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